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Abstract 



Searches for scalar top and bottom quarks have been performed with data collected by the ALEPH 
detector at LEP. The data sample consists of 21.7 pb~^ taken at ^/s = 161, 170, and 172 GeV and 
5.7 pb~^ taken at ^/s = 130 and 136 GeV. No evidence for scalar top quarks or scalar bottom quarks 
was found in the channels t — > cx, t — > hiv, and b — > bx- For the channel t ^ cx a limit of 
67 GeV/c^ has been set on the scalar top quark mass, independent of the mixing angle between the 
supersymmetric partners of the left and right-handed states of the top quark. This limit assumes a 
mass difference between the t and the x of at least 10 GeV/c^. For the channel t hiu the mixing- 
angle independent scalar top limit is 70 GeV/c^, assuming a mass difference between the t and the D 
of at least 10 GeV/c^. For the channel b — > bx, a limit of 73 GeV/c^ has been set on the mass of the 
supersymmetric partner of the left-handed state of the bottom quark. This limit is valid if the mass 
difference between the b and the x is at least 10 GeV/c^. 
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1 Introduction 



In the Minimal Supersymmetric Extension of the Standard Model (MSSM) |^ each Standard 
Model fermion has two scalar supersymmetric partners, one for each chirality state. The scalar- 
tops (stops) tR and tL are the supersymmetric partners of the top quark. These two fields are 
weak interaction eigenstates which mix to form the mass eigenstates. The stop mass matrix is 
given by ||: 



where mi^ and m^^ are the t^ and tL mass terms, at is related to the soft SUSY-breaking 
parameter Athj at = At — fi/ tan/? (where fi is the supersymmetric mass term which mixes the 
two Higgs superfields and tan (3 is the ratio between their vacuum expectation values) and rrit 
is the top quark mass. Since the off-diagonal terms of this matrix are proportional to rrit, the 
mixing between the weak interaction eigenstates may be large and the lighter stop could be the 
lightest supersymmetric charged particle. The stop mass eigenstates are obtained by a unitary 
transformation of the tR and tL fields, parametrised by the mixing angle 6i. The lighter stop is 
given by t = tL cos^^ + tR sin6't, while the heavier stop is the orthogonal combination. 

The stop could be produced at LEP in pairs, e"'"e~ — tt, via s-channel exchange of a virtual 
photon or a Z. The production cross section [Q] depends on the stop charge for the coupling 
to the photon and on the weak mixing angle and the mixing angle 6i for the coupling to 
the Z. When 6i is about 56° the lightest stop decouples from the Z and its cross section is 
almost minimal. At ^/s = 172 GeV, the maximum cross section is of order 1 pb for a t mass of 
60 GeV/ and is reached for 6i = 0°. 

The searches for stops described here assume that all supersymmetric particles except the 
lightest neutralino x and (possibly) the sneutrino D are heavier than the stop. The conservation 
of R-parity is also assumed; this implies that the Lightest Supersymmetric Particle (LSP) is 
stable. Under these assumptions, the two dominant decay channels are t — > cx or t hiu 0. 
The corresponding diagrams are shown in Figures ^ and Pd. The first decay can only proceed 
via loops and thus has a very small width, of the order of 1-0.01 eV 

The t hiu channel proceeds via a virtual chargino exchange and has a width of the order 
of 0.1-10 keV 0, where the largest width is reached for a chargino mass close to the stop mass. 
This decay dominates when it is kinematically allowed. Assuming equal mass sneutrinos u^, 

and the lepton flavour for this decay is determined by the chargino composition. If the 
chargino is the supersymmetric partner of the W the decays t — >• bez/g, t b/iz/^ and t hru^- 
occur with equal branching fractions. If the chargino is the supersymmetric partner of the 
charged Higgs the branching fraction of the decay t —>■ hrur is enhanced. In all of these cases, 
if the neutralino is the LSP the sneutrino can decay into {x^) but this invisible decay does not 
change the experimental topology. 

A possible third stop decay channel is the four-body decay t hfif2X- One such four-body 
decay of the t is shown in Figure |l|c. The rates of four-body decays are expected to be much 
smaller than that of the decay t ^ cx- 

The phenomenology of the scalar bottom (sbottom), the supersymmetric partner of the 
bottom quark, is similar to the phenomenology of the stop. In contrast to stops, sbottom mixing 
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Figure 1: Stop decay diagrams, (a) t cx- (b) t hiu. (c) t bfif2X- Decay (c) is not 
considered in this paper. 

is expected to be large for large values of tan because of the relation af, = Af, — fi tan p. When 
the sbottom mixing angle 6^ is about 68° the lightest sbottom decouples from the Z. Assuming 
that the b is lighter than all supersymmetric particles except the the b will decay as b — > bx- 
Compared to the t decays, the b decay has a large width of the order of 10-100 MeV. 

Direct searches for stops and sbottoms are performed for the stop decay channels t — > cx 
and t — > hiu and for the sbottom decay channel b —>■ bx- The results of these searches supersede 
the ALEPH results reported earlier for data collected at energies up to ^/s = 136 GeV 0]. The 
DO experiment [^] has reported a lower limit on the stop mass of 85 GeV/c^ for the decay into 
cx and for a mass difference between the t and the x larger than about 40 GeV/ c^. Searches 
for t — > cx, t — > hii> and b — > bx using data collected at LEP at energies up to y/s = 172 GeV 
have been performed by OPAL 



2 The ALEPH detector 

A detailed description of the ALEPH detector can be found in Ref. 0, and an account of 
its performance as well as a description of the standard analysis algorithms can be found in 
Ref. p]. Only a brief overview is given here. 

Charged particles are detected in a magnetic spectrometer consisting of a silicon vertex 
detector (VDET), a drift chamber (ITC) and a time projection chamber (TPC), all immersed 
in a 1.5 T axial magnetic field provided by a superconducting solenoidal coil. Between the TPC 
and the coil, a highly granular electromagnetic calorimeter (ECAL) is used to identify electrons 
and photons and to measure their energy. Surrounding the ECAL is the return yoke for the 
magnet, which is instrumented with streamer tubes to form the hadron calorimeter (HCAL). 
Two layers of external streamer tubes are used together with the HCAL to identify muons. 

The region near the beam line is covered by two luminosity calorimeters, the SICAL and the 
LCAL. The SICAL provides coverage from 34 to 63 mrad from the beamline while the LCAL 
provides coverage out to 160 mrad. The low angle coverage is completed by the HCAL, which 
occupies a position behind the LCAL and extends down to 106 mrad. The LCAL consists of 
two halves which fit together around the beamline; the area where the two halves come together 
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is a region of reduced sensitivity. This "vertical crack" accounts for only 0.05% of the total 
solid angle coverage of the ALEPH detector. 

The information obtained from the tracking system is combined with the information 
obtained from the calorimeters to form a list of "energy flow particles" |^. These objects 
are used to calculate the variables that are used in the analyses described in Section 4. 



3 Monte Carlo Simulation 

In the simulation of a stop signal, the most significant issues to be addressed are the treatment 
of the stop perturbative gluon radiation, hadronisation and decay. 

Since the stop is a scalar particle, the spectrum of gluon emission differs from that of a 
quark. The standard shower evolution programs would therefore need modifications to include 
the gluon emission from a spin- zero particle. However, as pointed out in Ref . , the difference 
between the average energy loss due to perturbative gluon emission off a spin-0 and a spin- 1/2 
particle is small ( ^ 10^^) and can safely be neglected within the approximations used by most 
shower Monte Carlo codes. 

The stop lifetime is longer than the typical hadronisation time of (9(10~^^ s), which 
corresponds to a width of (9(0.1 GeV). Stops therefore hadronise into colourless (tg) or (Iqq) 
bound states before decaying. This is incorporated in the generator by letting stops hadronise as 
if they were ordinary quarks according to the LUND string fragmentation scheme implemented 
in JETSET 7.4 [|T^. A Peterson fragmentation function |1T| is used to describe the stop 



fragmentation. The ei parameter in the function is scaled from b quarks following the relation 
ei = eb/n^/m? [O, with eb = 0.0035 and mb = 5 GeV/c^. Stop hadrons then decay 



according to a spectator model. The effective spectator quark mass M^s, which takes into 
account non-perturbative effects, is set to 0.5 GeV/ c^. The decay quark, c or b depending 
on the decay channel, is allowed to develop a parton shower to take into account hard gluon 
emission. At the end of the parton shower, a string is stretched among all coloured particles. 

A similar procedure is followed for the sbottom generator, taking into account the fact that 
the b lifetime is much shorter than the t lifetime. Depending on the b and x mass difference 
and coupling, the b can decay either before or after hadronisation. Two sets of b signal samples, 
one for each of these possibilities, were generated over the same range of mass differences. 

Signal samples were generated at ^/s = 130, 136, 161, and 172 GeV for various {mi,m^), 
{m^, m^) or (m^, rriy) masses. In these generations the mixing angle 6i or 6^ was set to zero; the 
selection efficiency depends on the value of the mixing angle, since changing its value changes 
the spectrum of initial state radiation. Two sets of t b£z> samples have been produced. The 
first set assumes equal branching fractions for the t decay to e, /x or r, while the second set 
assumes a branching fraction of 100% for the decay to r. All of these samples were processed 
though the full ALEPH detector simulation. 

The dependence of the selection efficiencies on the fragmentation parameters and on the 
mixing angle is discussed in Section 5. The effect of the short b lifetime on the b selection 
efficiency is also discussed in Section 5. 

Monte Carlo samples corresponding to integrated luminosities at least 100 times that of 
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the data have been fully simulated for the annihilation processes e^e~ ff and the various 
processes leading to four-fermion final states (e+e~ WW, e+e~ — > Weu, e"'"e~ Zee and 
e"'"e~ — * Z7*). The two-photon processes 77 £^i~ were simulated with an integrated 
luminosity about 20 times that of the data, while the two-photon processes 77 qq were 
simulated with an integrated luminosity about three times that of the data. 



4 Analysis 

Data collected at ^/s = 130, 136, 161, 170, and 172 GeV have been analysed, corresponding 
to integrated luminosities of 2.9, 2.9, 11.1, 1.1, and 9.5 pb~^, respectively. To account for the 
dependence on ^/s, all cuts are performed in terms of variables normalised to the beam energy. 

Two analyses are used to search for t production. The first one is sensitive to the decay 
t ^ cx while the second one is sensitive to the decay t —>■ hiu. Both channels are characterised 
by missing momentum and energy. The experimental topology depends largely on Am, the 
mass difference between the t and the x or ^- When Am is large, there is a substantial amount 
of energy available for the visible system and the signal events tend to look like WW, Wez/, 
Z7*, and gg(7) events. These processes are characterised by high multiplicity and high visible 
mass Mvis. When Am is small, the energy available for the visible system is small and the 
signal events are therefore similar to 77 qq events. The process 77 qq is characterised by 
low multiplicity, low Mvis, low total transverse momentum pt and the presence of energy near 
the beam axis. In order to cope with the different signal topologies and background situations, 
each analysis employs a low Am selection and a high Am selection. 

The values of the analysis cuts are set in an unbiased way following the iVgs procedure |T3 



In this procedure, the cut values are varied and applied to the background samples and the 
signal samples in order to calculate CT95, the expected 95% Confidence Level (C.L.) limit on 
the signal cross section. The final cut values used in the analyses are the ones which minimise 
cTgs. Cuts used to eliminate background from 77 — > qq events are not varied. Such events are 
difficult to simulate when they go into the low angle region of the detector. Conservatively, the 
values of the cuts used against 77 —>■ qq events are tighter than the values given by the iVgs 
procedure. 

The experimental topology of the process e^e~ ^ bb (b — » hx) is quite similar to that of 
the process e^e" — > tt (t — cx)- A common selection is therefore used to search for these two 
processes. 



4.1 Search for t ^ cx and b ^ bx 

The processes e+e~ — tt (t — cx) and e+e^ ^ bb (b ^ hx) are characterised by two 
acoplanar jets and missing mass and energy. Two selections are employed, one for the small 
Am case (Am < 10 GeV/ c^) and one for the large Am case (Am > 10 GeV/ c^). A common 
preselection is used against 77 — > qq events in both the low and high Am analyses. The number 
of charged particle tracks Nch must be at least four, Mvis must be larger than 4 GeV/ and pt 
(Figure ^) must be larger than 2%^/s, or A%^/s if the missing momentum points to within 15° 
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Figure 2: (a) pt for 77 — > and t — > cx at ^ — 161 GeV. The solid histogram gives the 
77 ~^ distribution, the dashed histogram gives the signal distribution for mi — 65 GeV/c^ 
and Am = 5 GcV/ c^, the dotted histogram gives the signal distribution for mi = 65 GeV/ 
and Am =15 GcV/ c^. The cut pt > 1%\fs is indicated by the arrow, (b) i?iso/-E'iepton for 99(7) 
and t — > \)iv at \/i =161 GeV. The solid histogram gives the 55(7) distribution, the dashed 
histogram gives the signal distribution for mi — 60 GeV / and Am = 20 GeV/c^. The cut 
-E'iso/-E'iepton < 4 is indicated by the arrow. In (a), the cut £'12° = has been applied. In (b), at 
least one identified electron or muon is required. Normalization for the plots is arbitrary. 



in azimuth from the vertical crack in LCAL. The polar angle of the missing momentum vector, 
Op^.^, must be greater than 18° and the energy detected within 12° of the beam axis, £^12°, 
must be less than 5%-v^. Both the acoplanarity and the transverse acoplanarity must be less 
than 175°. The acoplanarity is defined to be 180° for a back-to-back topology and is calculated 
from the momenta directions of the two event hemispheres, defined by a plane perpendicular to 
the thrust axis. The transverse acoplanarity is obtained by projecting the event onto a plane 
perpendicular to the beam axis, then calculating the two-dimensional thrust axis and dividing 
the event into two hemispheres by a plane perpendicular to that thrust axis. Both of these cuts 
are also effective against ^'^(7) background. 

4.1.1 Low Am selection 

Most of the cuts in the low Am analysis are designed to eliminate the remaining background 
from 77 qq events. The pt cut is reinforced by calculating pt excluding the neutral hadrons 
found by the energy flow algorithm and requiring it to be greater than 2%^/s. The pt is also 
calculated with only the charged particle tracks and required to be greater then l%^/s. These 
cuts eliminate 77 events that have a large pt due to spurious calorimeter objects; these objects 
can occur when soft tracks are not correctly associated with deposits in the ECAL or HCAL. 
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Such events are also eliminated by asking that the most energetic neutral hadronic deposit be 
less than 30% of the total visible energy £^vis- To eliminate 77 events that pass the pt cuts, 
£'12° must be equal to zero, Op^.^^ must be greater than 37°, ^thrust, the polar angle of the thrust 
axis, must be greater than 41°, and the missing mass Mmiss divided by Af^js must be less than 
25. Also of use is the fact that the missing momentum in 77 qq and qq{'~f) events can arise 
from neutrinos produced in semileptonic decays. When these decays occur within a jet, the 
resulting missing pt is not isolated. Signal events are therefore selected by requiring the energy 
in a 30° azimuthal wedge around the direction of missing pt to be less than 25%y/s. 

If a scattered electron from ^ qq process goes into an insensitive region of the detector, 
only a small fraction of the electron energy may be recorded. The missing electron energy can 
lead to a large measured pt, faking a signal. These fake signals can be eliminated by calculating 
the scattered electron angle ^scat from the pt, assuming the other electron to be undefiected, 
and by computing the angle ^point between the calculated electron direction and the closest 
energy deposit. The fake signals surviving the pt cut usually have a large value of ^scat; because 
the Pt imbalance is large, and a small value of ^point) because the calculated electron direction 
points to the energy deposit from the scattered electron. Both ^scat and ^point are incorporated 
into the analysis through the cut ^point > 60° — 10 x ^^scat- 

Additional cuts are used against the 77 t^t^ background. Most of the 77 — ^ t~^t^ events 
that survive the above cuts have four charged particle tracks from the decays r — > one-prong, 
T — > three-prong, and the low visible mass and high value of acoplanarity characteristic of 77 
events in general. In order to eliminate these events, any four-track event must have transverse 
acoplanarity less than 150° or visible mass greater than 20 GeV/ c^. As an additional safeguard, 
all four-track events are required to have a visible mass larger than 8 GeV/ regardless of the 
value of the transverse acoplanarity. 

The low Am analysis is completed by applying cuts against low mass WW, Z7*, and Weu 
events. A cut of thrust < 0.97 is effective against Z7* (with Z vv), while WW and Wei/ 
events are eliminated by requiring that E'vis be less than 26%-v/i. Events from the process 

WW — ^ iviTVr-, where the r subsequently undergoes a three-prong decay, are eliminated 
by requiring that the event mass excluding identified electrons and muons be greater than 
3 GeV/c2. 

4.1.2 High Am selection 

The main background in the high Am case comes from WW, Wez/, Z7*, and qq{'j)- Events from 
77 processes may still contribute to the background because they have a very large cross section 
and because detector effects may lead to extreme values for variables such as pt. Background 
from 77 is reduced by requiring that A'ch be greater than six and that pt be greater than 
5%^/s, or 7.5%^/s if the missing momentum points to within 15° of the vertical LCAL crack. 
Additional 77 events are removed by requiring that pt be greater than 20%£^vis- As in the 
low Am selection, it is necessary to guard against 77 events that have a large pt due to a 
missed association between soft tracks and calorimetry deposits. This is done by demanding 
that the total energy from neutral hadrons be less than 30%£'vis; this is relaxed to 45%ii^vis if 
the Pt calculated without neutral hadrons is greater than 3%-v^. Other cuts which are effective 
against 77 events are ^point > 5°, < 7.5%a/s and the total energy more than 30° away from 
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the beam greater than 30%-Evis- 

Finally, cuts against WW, Wei/, and Z7* are applied. Events from Z7* are eliminated by 
requiring that the thrust be less than 0.935. To eliminate WW events in which one of the W's 
decays leptonically, any identified electron or muon must have an energy less than 20% y^. In 
order to further reduce background from WW and Wei/ events, an upper cut is applied on the 
visible mass. The optimal value of this cut as determined by the Ng^ procedure depends on 
the mass difference of the signal sample considered. A hypothesis of Am = 15 GeV/c^ gives 
an optimal value of 0.315y^ for the M^is cut while a hypothesis of Am > 35 GeV/ gives an 
optimal value of 0.375^5 for the M^ig cut. 

The high Am selection changes as a function of Am through the Mvis cut. When this 
selection is applied to the data, the loosest M^^s cut is used. In the case that limits must be 
set, a candidate is counted for a given value of Am only if it has a visible mass less than the 
Mvis cut used for that value of Am. 

4.1.3 Selection efficiency and background 

To combine the low and high Am selections, three possibilities are considered: the low Am 
selection may be used, the high Am selection may be used, or both selections may be used. 
According to the Ng^ procedure the two selections should not be used simultaneously for any 
value of Am. For Am < 10 GeV/ c^, the low Am selection is used, while for Am > 10 GeV/ c^, 
the high Am selection is used. The t efficiencies are shown in Figure ^ while the b efficiencies 
are shown in Figure These b efficiencies are evaluated assuming that the b hadronises 
before it decays. 

For the low Am selection, the requirement that £"12° = results in an inefficiency due to 
beam-related and detector background. The size of this effect (~ 4%) has been measured using 
events triggered at random beam crossings and the low Am selection efficiency is decreased 
accordingly. 

The background to the low Am selection is dominated by 77 — » qq and 77 t^t^ and 
has a total expectation of 0.9 events (~ 40 fb) at ^/s = 161-172 GeV and 0.2 events (~ 30 fb) 
at a/s = 130-136 GeV. For the high Am selection, the background is dominated by WW, 
Wez/, Z7*, and qq{-f) at ^ = 161-172 GeV and by qq{-f) at ^ = 130-136 GeV. The total 
background expectation for this selection is 1.0 event (~ 50 fb) at a/s = 161-172 GeV and 0.2 
events (~ 30 fb) at y/s = 130-136 GeV, using the loosest value of the M^is cut. 

4.2 Search for t h£i) 

The experimental signature for t hiu is two acoplanar jets plus two leptons with missing 
momentum. The leptons tend to have low momenta, especially for low Am signals; when Am is 
8 GeV/ c^, the most energetic lepton often has a momentum between 1 and 2 GeV/ c. In order 
to identify electrons and muons, loose identification criteria based on the pattern of deposits 
in the ECAL and the HCAL have been applied. These loose criteria allow 1 GeV/ c electrons 
and 1.5 GeV/ c muons to be identified. Since low-momenta lepton candidates are often mis- 
identified pions, other analysis cuts must be used to keep the background at a low level. Two 
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Figure 3: Efficiencies as a function of Am. (a) Efficiency for a 65 GeV/c^ stop decaying as 

t — cx (soUd curve), a 50 GeV/c^ stop decaying as t — > cx (dashed curve) and a 60 GeV/c^ 
stop decaying as t — hiu (dotted curve), (b) Efficiency for a 60 GeV/c^ sbottom (solid curve) 
and a 50 GeV/ sbottom (dashed curve) decaying as b — > bx. 

selections are used, one for the small Am case (Am < 10 GeV/ c^) and the other for the large 
Am case (Am > 10 GeV/ c^). A preselection common to both the low and high Am selections 
is used against the 77 — qq background. It is required that Net be greater than six and M^^^ 
be greater than 8%y/s. It is also required that pt be greater than 1.25%y^, £'12° be smaller 
than 2 GeV, and ^point be greater than 50° — 20 x ^scat- In order to eliminate the radiative 7/7 
events in which a return to the Z peak has occurred, events with a longitudinal momentum 
greater than are rejected. 

4.2.1 Low Am selection 

If Am is small the visible energy is also small and both the jets and leptons are very soft. Since 
very soft leptons might not be identified, events with no electrons or muons are accepted. The 
main background arises from 77 — > qq. It is therefore required that £'120 = and that both 
^-Pmiss ^thrust be greater than 37°. An acoplanarity between 100° and 179° is also required. 
There must be at least one electron or muon with momentum greater than l%^/s, otherwise 
both the Pt cut and the two-dimensional cut in the 6'point-^scat plane are tightened: pt > 2%y^, 

^point > 115° — 20 X ^^scat- 

The WW background is eliminated by requiring that the missing mass be greater than 
S2.5%a/s and that the hadronic mass be smaller than 5%y/s if at least one electron or muon is 
identified. The qq{'y) events are rejected by requiring that the thrust be smaller than 0.9. 
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4.2.2 High Am selection 

For large mass differences at least one electron or muon with momentum between 2 and 
35 GeV/c is required. It is further required that E^so, the energy in a 30° cone around the 
direction of the electron or muon momentum (Figure ^d), be smaller than four times the electron 
or muon energy. If a second electron or muon is identified, i^iso is required to be smaller than 10 
times the electron or muon energy. If only one electron or muon is found, a tau jet is selected 
using the JADE algorithm with ?/cut = 0.001. This candidate r jet must have an energy smaller 
than 30 GeV, have less than 2 GeV of energy carried by neutral hadrons, and have an angle 
of at least 20° with the nearest jet. Finally, the missing mass is required to be greater than 
25%v/i. 

To reinforce the 77 —>■ qq rejection further cuts are needed. It is required that Op^_^_^^^^ be 
greater than 18°, that the transverse acoplanarity be smaller than 176° and that the acollinearity 
be smaller than 174°. If only one electron or muon is identified the hadronic neutral mass must 
be smaller than 30%£'vis and the cuts on Op^.^^ and pt are tightened: Op^^.^^ > 26°, pt > 3%y/s. 

The WW background events are eliminated by requiring that M^is be smaller than 74% 
and that the hadronic mass be less than 37%^/s. It is also required that the quadratic mean 
of the two inverse hemisphere boosts {^J {{mi/ EiY + {jn^j E-if^ jl with m\^2 and E\^2 the two 
hemisphere masses and energies) be greater than 0.2. The remaining gg(7) background is 
reduced by requiring that the thrust be smaller than 0.925. 

4.2.3 Selection efficiency and background 

The low and high Am selections are combined using the same procedure as in Section 4.1.3. In 
contrast to the situation for the t — *• cx channel, the smallest value of CT95 is obtained when the 
low and high Am selections are used simultaneously. This is true for all values of Am. Shown 
in Figure ^ is the efficiency assuming equal branching fractions for the t decay to e, [i or r. 
If the branching ratio to r is 100%, the efficiency is about 35% for a Am between 10 and 35 
GeV/ (? . As is the case for the t — >■ cx channel, the inefficiency caused by the beam- related 
and detector background is taken into account. 

Most of the background comes from the high Am selection and is dominated by 55(7) at 
v/i = 130-161 GeV and by WW and gg(7) at ^/s = 170-172 GeV. A total of 0.8 events (~ 30 fb 
at 161 GeV and ~ 50 fb at 172 GeV) are expected qX = 161-172 GeV while 0.2 events 
(~ 30 fb) are expected at ^/s = 130-136 GeV. 

5 Systematic Uncertainties 

The systematic uncertainty on the t and b selection efficiencies comes mainly from the limited 
knowledge of t and b physics (hadronisation and decay). Uncertainties related to detector 
effects, to the size of the signal samples, and to the parameterisation of the signal efficiencies 
are also considered, and for the t — > hii) analysis the effects of lepton identification are taken 
into account. The physics model used in the generators is described in Section 3; the systematic 
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effects are studied by varying the parameters of the model and checking the resultant effect on 
the efficiency. 

The change in the efficiency due to the systematic effects is shown in Table 3. When Am 
is small, the uncertainties associated with the t and b physics are relevant. The largest change 
in the low Am efficiency comes from the variation in M^s- This variation changes the invariant 
mass available for the hadronic system and thus the multiplicity and event shape. To quantify 
these effects, Meg is varied from 0.3 GeV/ to 1.0 GeV/ c^, a range much larger than that 
implied by low energy measurements. When Am is large, the selection efficiencies are insensitive 
to the values of the parameters, changing by only ~ 2% relative even for Mcs = 2 GeV/ c^. 

The fragmentation parameters are varied over a range suggested by LEPl measurements. 
In the case of the error is propagated from eb according to the formula described in Section 3, 
and for the t — * hii) channel eb is varied simultaneously with e^. Similarly, for the b bx 
channel eb is varied simultaneously with ej^. For the large Am case the fragmentation parameters 
are varied more drastically, but even drastic variations have little effect on the efficiency; for 
instance, when ei = eb, the relative change in large Am t efficiencies is only ~ 2%. 

The systematic effect of varying the mixing angles is quantified by evaluating the efficiencies 
on a set of t samples generated with 6i = 56° and on a set of b samples generated with 6^ = 
68°. For these values of mixing, the stops and sbottoms decouple from the Z and the change 
in efficiencies due to differing amounts of initial state radiation is maximal. 

The structure of the matrix element in the semileptonic decay t hiu is also considered. 
Two sets of t — s> hiu signal samples are generated. One set includes the the matrix element, 
treated as in Reference 2, while the other set employs a phase space decay model. Including 
the matrix element increases the efficiency of the t hiu selection by about 5% relative with 
respect to the phase space decay model. Conservatively, the phase space decay model is used 
to obtain the t hiu efficiencies. 

The effect of the relatively short b lifetime has been checked by comparing the two sets of b 
signal samples. Higher efficiencies are always obtained from the set in which the b decays before 
hadronisation. The lower efficiencies, obtained under the assumption that the b hadronises 
before decay, are taken as the actual efficiencies; this helps ensure that any limits set on b 
production will be conservative. 

The size of the signal samples, 1000 events, leads to a relative uncertainty of less than 2%, 
while the parameterisation of the signal efficiencies leads to an additional relative uncertainty 
of ~ 2%. The total statistical uncertainty associated with the Monte Carlo signal simulation 
is therefore ~ 3% relative. 

Detector effects have been studied for the variables used in the analyses. Events in the 
data from gg(7) final states are selected with a loose set of cuts and compared with the qqi^j) 
Monte Carlo. All of the relevant variables, such as pt and ^point, show good agreement. The 
lepton isolation and the lepton identification, which are crucial for the t — ^ hli) analysis, are 
also considered. The lepton isolation shows good agreement between gg(7) Monte Carlo and 
data, while the lepton identification is found to lead to a 3% systematic error. 

The systematic errors are incorporated into the final result using the method described in 
Reference flH . 



10 



Table 1: Summary of relative systematic uncertainties on the t and b selection efficiencies. 
The ranges of variation are those used for the low Am case. 



Systematic Uncertainties (%) 



Type 


t > 


cx 


D — > 


bx 


4- V 
t > 






High Am 


Low Am 


High 


Low 


High 


Low 


Meff (0.3-1.0 GeV/c2) 


3 


10 


4 


11 


3 


15 


ei, eb(eb0.002 - 0.006) 


2 


2 






2 


2 


eb,eb(eb0.002 - 0.006) 






1 


2 






ec (0.02-0.06) 


3 


7 










9i (0°-56°) 


1 


3 






2 


1 


9^ (0°-68°) 






3 


2 






Monte Carlo statistics 


3 


3 


3 


3 


3 


3 


detector effects 


negl. 


negl. 


negl. 


negl. 


3 


3 


TOTAL 


6 


13 


6 


12 


6 


16 



6 Results 

One event is selected by the t cx, b ^ bx selection, while no events are selected by the 
t hiu selection. The candidate event is selected at ^/s = 161 GeV; its kinematic properties 
suggest the process e"'"e~ — > Z7* — >• i/z/r+r" as a Standard Model interpretation. Since only 
a single event is selected, it is appropriate to set lower limits on the masses of the t and b. 
Figures ^a, ^d, and ^ give the 95% C.L. excluded regions for the channel t cx- For this 
channel, the ^^-independent lower limit on mi is 67 GeV/ c^, assuming a mass difference between 
the t and the x of at least 10 GeV/c^. Figures ||a, Qd, and |^c give excluded regions for the 
t hiu channel, assuming equal branching ratios for the t decay to e, /i, r. In this case, the 
^^-independent lower limit on m^ is 70 GeV/ c^, assuming a mass difference between the t and 
the z> of at least 10 GeV/c^. 

Figure ^d gives the excluded region in the (Am, mi) plane for the t — > hiV channel, assuming 
a branching ratio of 100% for the t decay to r. A 6't-independent lower limit of 64 GeV/ is 
set on mi in this case, again assuming a mass difference between the t and the u of at least 
10 GeV/c^. 

Figures |a, §3 and ^ give the excluded regions for the b decay b hx- A lower limit of 
73 GeV/c^ is set on m^ , assuming that 9^ is 0° and that the mass difference between the b 
and the x is at least 10 GeV/c^. Figure ||b shows that 6'j^-independent limits are not set. 
When decoupling from the Z occurs, sbottoms can only be produced through photon exchange 
and the cross section for the b (charge —1/3) is four times lower than the cross section for the 
t (charge +2/3). 
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7 Conclusions 



Searches have been performed for scalar top quarks at = 130-172 GeV. A single candidate 

event, selected at a/s = 161 GeV, is observed in the t cx channel while no events are observed 
in the t h£u channel. This is consistent with the background expectations of 2.3 events for 
the t ^ cx channel and 1.0 events for the t hiu channel. 

A 95 % CL. limit of mi > 67 GeV/ is obtained for the t — cx channel, independent of 
the mixing angle and valid for a mass difference between the t and the % larger than 10 GeV/ c^. 
For the t — > h£u channel, the ^^-independent limit is mi > 70 GeV/c^ if the mass difference 
between the t and the i> is greater than 10 GeV/c^ and if the branching ratios are equal for 
the t decays to e, and r. 

A limit is also obtained for the b decaying as b —> hx- The limit is m^ > 73 GeV/c^ for 
the super symmetric partner of the left-handed state of the bottom quark if the mass difference 
between the b and the x is greater than 10 GeV/ c^. 
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Figure 4: Excluded regions assuming t cx- (a) Excluded region in the vs mi plane, 
including the region excluded by the DO collaboration, (b) Excluded region in the mi vs 9i 
plane, (c) Excluded region in the mi vs Am plane. In (a) and (c), excluded regions are given 
for 0°, corresponding to the maximum t-Z coupling, and for 56°, corresponding to the minimum 
t-Z coupling. 
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Figure 5: Excluded regions assuming t hiu. (a) Excluded region in the m/> vs mi plane, (b) 
Excluded region in the mi vs 6i plane, (c) Excluded region in the mi vs Am plane. In (a), 

(b) and (c) equal branching fractions for the t decay to e, or r are assumed, (d) Excluded 
region in the mi vs Am plane, assuming a branching ratio of 100% for the t decay to r. In (a), 

(c) , and (d), excluded regions are given for 0°, corresponding to the maximum t-Z coupling, 
and for 56°, corresponding to the minimum t-Z coupling. Also shown in (a), (c), and (d) is the 
excluded region from LEPl, obtained from the measurement of the Z lineshape. 
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Figure 6: Excluded regions assuming b —>■ hx- (a) Excluded region in the vs plane, (b) 
Excluded region in the mj^ vs 6^ plane, (c) Excluded region in the vs Am plane. In (a) 
and (c), excluded regions are given for 0°, corresponding to the maximum b-Z coupling, and 
for 40°. 
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